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1^ I Abstract. The analytical solution of the three-dimensional linear pendulum in a 
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^ . and expressed in terms of initial conditions. This result offers the possibility of treating 

C^ I Foucault and Bravais pendula as trajectories of the very same system of equations, 
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1. Introduction 

During the first year(s) of undergraduate studies students are constantly warned to work 
in an inertial frame, as only in this case Newton's second Law can be safely applied. 
Since it rotates around its axes (even ignoring its motion around the sun), a reference 
frame fixed to the Earth is definitively non inertial. However, for most mechanics- related 
results this fact can be overlooked, as this rotation is not fast. But there are certain 
circumstances when non inertial reference frames, specially rotating frames, look like an 
interesting thing to use. This is specially so, when trying to study movement on Earth 
surface or when looking to prove or measure Earth motion itself through a mechanical 
experiment. 

The first proof of Earth's rotation, and perhaps the most simple and elegant, 
is indeed the Foucault pendulum experiment [I] performed in 1851. A handwriting 
annotation by Foucault founded among his papers reads: Mercredi, 8 Janvier, 2 heures 
du matin: le pendule a tourne dans le sens du mouvement diurne de la sphere celeste [2]. 
Namely, the experiment succeeded on January 8th at 2 hours a.m. It was carried out in 
his home's cellar, d'Assas street in Paris, using a 2 m. long pendulum [3]. Afterwards, 
in February, the experiment was repeated at a larger scale with a 11 m. long pendulum 
in the Paris Observatory [I] under the permission of Arago. It is worth noting the deep 
acknowledgement Foucault expresses toward his assistant Froment. 

Nowadays it is widely known the public exhibition of the pendulum held in 
March, on the request and backing of the prince-president of the Republic. Gracefully 
suspended from the ceiling of the dome of the Pantheon in Paris, it consisted of a 28 
kg brass-coated lead bob attached to a 67 meter long steel aircraft cable. The plane of 
the pendulum's swing appeared to rotate clockwise 264 degrees every day (11 degrees 
per hour), sweeping a circle of 16 m. diameter. Since no rotational forces act on the 
pendulum, it is obvious that it is the Earth beneath that is actually rotating. However 
sizeable effects, at naked eye, take several hours to become apparent. We will go back 
to this point later. 

On the very same year (and amazingly also in Paris) Bravais proposed an alternative 
way to observe the rotation of the Earth, using also pendula [H [5]. His approach was 
however different and has been review recently in this journal |6j. He explored the 
difference in the period of oscillations of right and left handed pendula. This fact, as 
we are going to show later, opens the door to observe the rotation of the Earth in a few 
minutes. 

Yet, despite not receiving a unified treatment, Foucault and Bravais pendula share 
the same equations of motion, differing only through initial conditions. Therefore they 
may be considered, in all senses, as two particular cases of a large variety of oscillations 
modes of a pendulum in a non inertial reference frame. 

Foucault type oscillations are obtained when the bob is moved away from the 
equilibrium position and then released with zero initial velocity. In an inertial reference 
frame the pendulum would oscillate in a fixed plane. Due to the Earth rotation, we 
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observe this plane to rotate. Thus the trajectory followed by the bob produces the 
well-known beautiful patterns, which are the proof of Earth's rotation. Instead, Bravais 
pendulum refers to conical oscillations of the bob. To this end, the bob is released 
from a non equilibrium point with a precise tangential speed in such a way that the 
projection of the trajectory on the horizontal plane describes simply a circumference. 
What renders Bravais pendulum interesting is that the trajectories are not invariant 
under sign reversal in the initial tangential velocities, as a consequence of Earth's 
rotation. It is just this phenomenon that led Bravais to propose the pendulum as 
an experimental demonstration of Earth's rotation. 

Being just two particular sets of initial condition for the same system of equations, 
it becomes obvious that there is a whole plethora of systems which presents interesting 
features worth looking at. Specially knowing that extremely precise initial conditions 
can be quite challenging to achieve from an experimental point of view. 

One of the purposes of this paper is to give an explicit solution of the linear 
pendulum in terms of initial conditions. We find it of particular interest for students 
whose training in solving linear systems of differential equations is concurrent with their 
course in mechanics. 

The structure of the paper is the following. In Section [2] we collect the equations 
of motion. In Section |3] we review some common approximate solution presented in 
textbooks. Next, in Section H] the analytical exact solution in three dimensions is 
obtained. Section |5] connects our results with a nowadays almost forgotten theorem 
by Chevilliet. In Section [6] a comparison two- and three-dimensional treatment of the 
Foucault pendulum is carried out. A previously unnoticed pattern in the attractor is 
reported. Eventually, Section [7] contains a discussion of our results. 

2. Equations of motion of the pendulum for small amplitudes 

We will follow the careful exposition in [7] to obtain the equations of motion of a 
pendulum with small oscillations, referred to a rotating reference frame. The starting 
point is the Newton equation for the bob of the pendulum 

mf=T — mgk, (1) 

where the dot stands for time derivative, r is the position vector of the bob of mass 
m with respect to an inertial system of reference, and T is the tension of the wire. In 
this instance, k is the unit vector directed vertically upward, which follows the radial 
direction of the Earth, assumed to have spherical symmetry. This remark will be of 
interest momentarily. Eventually, g is the corresponding intensity of gravity. 

The analysis of Foucault pendulum is carried out in a rotating reference system 
fixed to the Earth surface. We will consider that the Earth turns with constant angular 
velocity 

n = -incos/3 + knsmf3, (2) 
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Figure 1. 

pendulum. 



Reference system attached to the rotating Earth used to study the 




Figure 2. For a simple pendulum with small amplitude approximations gives 
Ti — — T i/l, with i — x,y,l ~ z and T ~ rag. 



in terms of the latitude /3, and Vt = 27r/86400 rd/s = 0.000073 rd/s. The unit vector i is 
tangent to the local meridian and points toward the equator, in the northern hemisphere. 
The third unit vector, j, is directed eastward (see Figure [1]). 

In a rotating reference system with origin located at the center of the Earth, time 
derivatives have to be replaced with more involved expressions. In particular, if fi = 0, 
the acceleration r becomes f+2Qxf+Qx [Q x f). The term 2Q x f is the Coriolis 
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acceleration, and i7 x (fi x r) stands for the centrifugal acceleration. Notice that both 
contributions are first and second order in ^2, respectively; and since fi <^ 1, they can 
be considered small corrections. 

If we move to a frame of reference fixed on the surface of the Earth at position R, 
then we have to replace r with r + R, in the transformation above. The equation of 
motion in the rotating system reads then 

■mr = T — mgk — m^ x {^ x R) — 2m^ x r — mVl x {Vl x r), (3) 

where, once again, r is now referred to the frame of reference on the surface of the 
Earth. The term g^ = gk + niQ x [Q x R) defines the so-called effective gravity, which 
takes into account the slight effect of centrifugal acceleration (which is order f2^). The 
vector ge determines the local plumb line and allows us to simplify the structure of 
([3]) by redefining the unit vector k so as to follow the direction of ge- This way, g^ 
has only one non-vanishing component (its value varies slightly with latitude) and is 
everywhere orthogonal to the surface of the Earth. We are implicitly stating that the 
Earth is not a sphere but its real shape is given by the equipotential surface resulting 
from the interplay between the gravitational and centrifugal interactions, the geoid. The 
definition of the unit vectors i,j does not change. They are, of course, orthogonal to k. 
These considerations simplify ([3]) to the form 

mr = T — mgk — 2mQ x f — mQ x [Q x r), (4) 

The linear approximation for the pendulum equations corresponds to assume that 
the tension vector is given by (see Figure [2]) 

f = -!^[xl + yj + {I ~ z)k], (5) 

valid for small amplitudes. Here, I is the length of the pendulum, and the origin of 
coordinates is at its supporting point. Introducing the notation u = ^/g/l, we get from 
(H and dS]) 

f =a;V-2fi X f-il X (^ X f), (6) 

or, explicitly in components 

X = — iJ^x + 2sfi?/ + s^fi^x + scVi^z, 

y = - uj^y - 2sVLx - 2dlz + fi^y, (7) 

z = — uPz + 2cVLy + csVl^x + (?Vl?z, 

where we have used the notation 

s = sin/3, c = cos/3. (8) 

We remind the reader that in these equations the z axis is vertical upward, the x axis 
is north-to-south, and the y axis is west-to-east. 

The approximation in ([5]) is crucial from a practical point of view, for it leads to 
the linear differential system ([7]), which is homogeneous with constant coefficients, and 
therefore solvable. One standard solution method consists in transforming ([7]) into a 
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linear differential equation of sixth order. Else, one can alternatively transform ([7]) into 
a system of first order and dimension six. Both procedures are indeed cumbersome and, 
to the best of our knowledge, no explicit exact solution is found in textbooks. 

3. Approximate trajectories in two dimensions 

In addition to the assumption of small amplitude oscillations implicit in ([5]) , two further 
approximations are usually carried out in order to obtain approximate trajectories. 
Firstly, one can assume Vt <^ u, and then ignore the centrifugal acceleration, namely 
the quadratic terms Vl? in ([7]). This provides a remarkable simplification. A further 
approximation considers that the pendulum length / is large enough so as the motion 
of the pendulum takes place in a plane. Thus z, i, z are dropped in ((Tj) and the 
equations of motion of the planar pendulum in the small amplitude approximation 
(without centrifugal acceleration) reduce to the two-dimensional system 

X = — u^x + 2sQy, 

y = - w^y - 2sQx. (9) 

Due to the symmetry of these equations a convenient way to solve them is to introduce 
the complex representation Z = x + iy. Then the system ([9]) may be written down in 
terms of a sole complex equation 

Z + 2iVLZ + u'^Z = 0, (10) 

with Cl = sQ. This is a linear homogeneous differential equation of second degree 
whose solution reads Z{t) = y4exp(A+t) + i?exp(A_t); and A,B are complex arbitrary 
constants. Here X± = —i{Q, ± ^),^ = y uj"^ + fi^, are the (complex) roots of the second 
degree equation A^ + 2i^\ + w^ = 0. The complex solution becomes 

Z{t) = exp{-int)[Aexp{i^t) + B exp{-i^t)]. (11) 

The four real quantities contained in the two complex constants A, B are determined 
so as x(t),y{t), satisfy the given initial conditions x{0),x{0),y{0),y{0). This requires 
some further algebra, but just to interpret the structure of the solutions let us think of 
initial conditions such that A and B are complex conjugate of each other. In that case 
the terms in square brackets stand for a real function and the solutions x{t),y{t) may 
be simply written down as 

x(t) = acos{Qt) cos{ijjt + (p), 

y(t) = a sm{Qt) cos {out + (p), (12) 

in terms of an amplitude a, a phase 0, and where, in addition, we have approximated 
C, ^ u. Thus the trajectories as a function of t are the product of a fast oscillation of 
frequency u describing the natural motion of a simple pendulum, and a slow oscillation 
of frequency Q which originates a rotation of the pendulum oscillation plane. 

By using elementary matrix methods, the general solution of iQ may be expressed 
in terms of initial conditions. To get rid of the complex representation in (ITTl) we write 
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r, u,j, 6j., &i- Next, using 



A ^ ttr + itti^ B ^ h,f. + hi, in terms of new four real parameters a. 
Euler formula in Z{t) and Z{t), equation fITT]) yields after some algebra 

I x{t) \ I cos at — sincrt cos pt sin pt 

x{t) —a sin at —p cos at — p sin pt p cos pt 

y{t) sin at cos at — sin pt cos pt 

V 2/(^) / \ crcoscrt —cr sincrt — p cos pt —psinpt j 

where cr = ^ — r2,p = .^ + i7, with ^ already defined above. Equation ( lT3ll defines the 
time dependent matrix M{t). At time t = we have 

/ a;(0) \ / flr \ 

(14) 
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i;(0) 
1/(0) 
V 2/(0) / 
Therefore, we can now give the trajectories in terms of initial conditions 

/ x(0) \ 

i;(0) 

2/(0) 
V 2/(0) / 

\ 
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The explicit solution for the planar, small amplitude pendulum, with Coriolis and 
without centrifugal accelerations reads then 
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y{t) 
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\ m 1 


straightforward computation yields 


/ p 1 


M-i(O) = 


-1 p 
a - 








1^ 1 a 



(16) 



x(t) = - [(.^ cos fit cos ,^t + r2sinntsin^t)xo + (cosr2tsin^t)i;o + 



(17) 



y(t) = [{C, sin Clt cos ^t + Cl cos Clt sin C,t)xQ 



;i8) 



(^ sin fit cos ^t — ficosfitsin,^t)yo + (sinfitsin,^t)yo], 

(sinfitsin^t)i;o + 
(^ cos fit cos ^t + fi sin fit sin $,t)yo + (cos fit sin ^t)yo] , 

where we use the notation a;(0) = xo,y{0) = yo,x{0) = xo,y{0) = ijQ. 

Equations ( TT7I) and ( TT8l) allow us to illustrate the shape of trajectories for different 
sets of initial conditions, and in particular those that give rise to the Foucault and 
Bravais pendula. The Java applet freely downloadable at [8] allows an easy visualization 
of the trajectories for arbitrary values of the parameters. 

Next we illustrate the pendulum trajectories with three different choices of initial 
conditions. We have chosen the following parameter values: Length of the wire L = 100 
m., which gives the angular velocity to = \/ gjL = 0.3162 rd/s; rotating reference system 
velocity fi = a;/10; and latitude /3 = 7r/4. The Earth rotation velocity is artificially 
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speeded up for the sake of clarity in the figures. In all the figures, the solid dot stands 
for the initial position. The arrow establishes the direction of turn. 

3.1. Foucault pendulum in 2D 

The Foucault pendulum corresponds to the simplest possible choice of initial conditions. 
It is obtained when we release the bob with vanishing initial velocity (in the Earth 
attached system) from, say, the point {xq 7^ 0,yo = 0). We illustrate in Figure El left 
panel, a trajectory in the plane xy. The initial values are: Xq = l,yo = 0, Xq = ijo = 0. 
Also, in the right panel of Figure [3] we plot the corresponding trajectory in the plane of 
velocities xy. 

This case plainly and elegantly combines the motion of a single, planar pendulum 
with the precession of its plane of oscillation. The deflexion of the plane of oscillation 
originates in the Coriolis acceleration. Thus, the term 2^2 x r deflects the trajectory to 
the right of the velocity vector in the northern hemisphere. Hence, every half oscillation 
of the pendulum produces an incremental effect. The drawback of this set up is that, as 
we stated before, the effects take several hours to become evident. For a Foucault 
pendulum located at the latitude of Paris for example, it takes about 32 hours to 
complete a precession cycle. In order to plot neat trajectory patterns we have used a 
much more greater angular velocity than the one corresponding to the value of Earth's 
rotation. 

3.2. Bravais pendulum in 2D 

The Bravais pendulum or conical pendulum, is doubtlessly the second simplest choice 
as regards to the initial conditions. It is obtained when we release the bob from a point 
(xo 7^ 0,?/o = 0) with an initial tangential velocity (xq = 0, yo = xoiuj — Vt)). Then the 
pendulum describes conical oscillations and the planar trajectory is a circumference. 

If we reverse the initial velocity, (xq = 0, yo = —xq{uj — i7)), the trajectory is only 
conical-like. The planar projection is not a circumference any longer, but rather an 
open trajectory. This is due to the effect of the rotation of the Earth on the motion of 
the pendulum bob. 

This is illustrated in Figure HI The black solid line corresponds to the conical 
pendulum released at t = with the initial conditions: xq = l,yo = 0,i;o = 0,yo = 
0.13142. Whereas the red dotted line describes the trajectory started with reversed 
initial velocity. In Figure HI the left panel corresponds to the trajectory in the plane 
and in the right panel we plot the corresponding trajectory in the plane of velocities xy. 

Due to Earth's rotation, the left handed and the right handed pendula take slightly 
different times to complete one turn, besides the slight difference in the shape of their 
paths. The delay is accumulative and their mutual phase difference becomes sizeable 
after completing a moderated number of turns. It therefore provides an easy and fast 
proof of Earth's rotation, at least in an ideal experiment with two simultaneous pendula. 
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Figure 3. Foucault pendulum in two dimensions. Trajectory in phase plane (left 
panel) and in the velocity phase plane (right panel). The solid dot stands for the 
initial position. For the numerical values of the parameters see main text. 
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Figure 4. Bravais pendulum in two dimensions. Trajectories in phase plane (left 
panel) and in the velocity phase plane (right panel). The black solid line represents 
the trajectory when the bob turns in the same direction as il. The red dotted line 
corresponds to a motion started with the reversed initial velocity. The solid dot stands 
for the initial position. For the numerical values of the parameters see main text. 



3.3. General pendulum in 2D 

It is also interesting to observe the trajectory when the initial velocity has two non- 
vanishing components. In Figure [5] (left panel) we plot the trajectory corresponding to 
initial conditions Xo = l,|/o = 0,Xo = — a;o(w — fisin/3),?/o = Xo(w — fisin/3). In the right 
panel of Figure |5] we represent the corresponding trajectory in the plane of velocities 
xy. 
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Figure 5. General pendulum in two dimensions. Trajectory in phase plane (left 
panel) and in the velocity phase plane (right panel). The solid dot stands for the 
initial position. For the numerical values of the parameters see main text. 



4. The linear pendulum in three dimensions: Exact solution 

Let us consider the equations of motion (j4]) of the pendulum in three dimensions without 
any further approximation, namely, we keep both Coriolis and centrifugal accelerations 



d T —* dv — * — * 

-^ + 2Qx \-Q X (Q xf) + uj'^r 

dt^ dt 



0. 



(19) 



The vector products above, being linear operations, admit a matrix representation. For 
instance 

/ -sin/3 \ 

sin /3 cos /3 

-cos/3 



VLx r = VL 



\ 



f = QAf. 



(20) 



/ 



Equation ( l20l) defines the matrix A of the cross product associated to the unit vector 
Cl. By the same token, we have Q x f= QAf, and Q x [Q x f) = Q'^A'^f. Now (fT9|) may 
be written down in matrix form 



dt'^ dt 



0, 



(21) 



which we solve next. 

We commence with a common exercise posed to undergraduate students in courses 
of ordinary differential equations. They are asked to transform the scalar second order 
differential equation 

pi + a{x)^ + b{x)y = 0, (22) 

into its canonical form, namely the one where the coefficient of the first order derivative 
vanishes. To this end one proposes a factorization of the dependent variable y{x) as 
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y{x) = u{x)v{x), with v{x) the new dependent variable and u{x) to be determined later 
on. Substitution in fl22|) yields 

(Pv f 2 du \ dv 1 f d'^u ^y h\ n (o'i\ 

dx^ \udx J dx u \ dx^ dx I 



If we chose 



u{x) = exp 



1 f 
-- / a{s)ds 



then we eventually get the canonical form 



d^ 
dx"^ 



b{x) - -{a{x)f - 2^"(^) 



(24) 



0. (25) 



-^ + u;V=0. (29) 



To put forward now the question of finding the canonical form of matrix equation 
(pTj) is in order. Following a similar procedure as above, we propose the substitution 

f(t) ^ f/(t)g(t), (26) 

in fl2T|) . where [/ is a square real matrix function of dimension three. Then, the canonical 
matrix differential equation, analogous to (!23|) . reads 

^U + nAU = 0, (27) 

dt 

Its solution reads 

U{t) = Cexp{-QAt), Cen^''^. (28) 

Since the particular value of the constant matrix C is irrelevant for our purposes, we 
choose the simplest option, namely C = I. As U{t) is not singular we get, after some 
algebra, the canonical form 

dt' 
Its solution in terms of initial conditions reads 

q{t) = cos{ujt) q{0) + - sin(wt) ^(0). (30) 

Thus, it turns out that the idea of finding the canonical form of the equations of 
motion has led to an uncoupled problem. Moreover, equation ( 129|) has an important 
dynamical interpretation. It tells us that U{t) transforms the dynamics into three free 
uncoupled oscillators. We know that this picture corresponds to an analysis in the 
inertial coordinate system fixed in space (i.e., not turning with the Earth). In this 
reference system the motion of the pendulum is observed as the composition of three 
free harmonic oscillators along the axes x, y, z. 

At this point a comment about the rhs of ( !28|) is in order. The exponential of a 
square matrix B is defined by the Taylor expansion: 

expB = I + B + ^B' + ^B'... (31) 

In general, the problem of finding the exponential of a matrix may be involved, but 
efficient methods do exist [21 [ID]- However, due to the symmetry of the matrix A, in 
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our case the computation may be easily carried out by inspection of the powers A". It 
is easy to verify that 

A^ = -A, A^ = -A^. (32) 

So, recurrent substitution in the Taylor expansion yields 

exp{nAt) = I - sm{nt)A + (1 - cos nt)A'^, (33) 

which provides us with an explicit matrix representation of U{t). 

To recover the solution in the rotating reference system fixed to the Earth we use 
equation ( l26l) 

f(t) = U{t)q{t) = e-^^* 



cos(a;t) g(0) H — sin(ci;t) g(0) 

CO 



(34) 



We still have to establish the relationship between the initial conditions in the inertial 
reference system: ^(0),g(0); and those in the rotating reference system fixed to the 
Earth: f(0),f(0). Equation ( 1M|) and its derivative at time t = give 

f(0) = g(0), f(0) = <f(0) + fiAg(0). (35) 

The exact solution in terms of initial conditions, in the system of reference that rotates 
with the Earth, after some algebra is found to be 

f(t) = -{[t^(l - cos fit) cos wt - fi sin cut sin fit] A^ r (0) + (1 - cosfit) sinwt A^ r(0) 

+ [fi cos fit sin ujt — OJ sin fit cos cut] A r(0) — sin fit sin ut A r(0) 

+ ci;coscutr(0) + sinci;tr(0)}. (36) 

Taking into account trigonometric identities for the product of sin and cos functions we 
can conclude that the motion is given by the superposition of three frequencies: cu ± fi, 
and u. This is at variance with the approximate trajectories in two dimensions ( 1T7|) and 
( 1T8|) where only the two frequencies: ^ ± fi ~ cu ± sfi, are present. In the case of the 
Earth rotation, it must not be relevant, for we should be able to resolve between two 
and three peaks in the power spectrum of, say, x{t)\ where fi ^ cu. 

Since that A is the matrix representation of the vector product associated to the 
unit vector fi, we can also write 

r{t) = — {[cu(l — cos fit) cos cut — fi sin cut sin fit] fi x (fi x r(0)) 

cu 

+ [(1 - COS fit) sin cut] fi X (fi X f(0)) 

+ [fi cos fit sin cut — cu sin fit cos cut] fi x r*(0) 

- [sin fit sin cut] fi x f(0) + cu cos cut f(0) + sincut f(0)}. (37) 

Notice that in the limit fi — ?■ 0, the equation above gives the expected result f{t) = 
cos cut r(0) + -sincut r(0). 

In terms of components the exact solution reads 

x{t) = I [c^ + cos (fi t) s^] cos (cu t) + — sin(fi t) sin(cu t) s^ I a;(0) 
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cos (fi t) sin (u t) — sin (fi t) cos (a; t) 



sy{0) 



+ 



u 
sinut 



sin(r2 1) sin(a; t) — [I — cos (fi t)] cos {u t) > sc z{0) 



UJ 



[c^ + COS (fi t) s^]x(0) + sin (fi t) s y{0) - [1 - cos (fi t)] sc i(0) [ (38) 



yit) 



n 



u 



COS [Q t) sin (w t) — sin [Vt t) cos {uj t) 



[sa;(0) + cz(0)] 



fi 



sin cut 



+ <! — sin(fi t) sin((X' t) + cos (f2 1) cos (w t) > y(0) 

sin {VL t) s x(0) - cos {Vt t) y(0) + sin {Vt t) c i(0) 



a; 
z{t) = \ — sin(f2 1) sin(a; t) — [1 — cos {Vl t)] cos (w t) \ sc x{0) 



(39) 



+ 



sin {Q t) cos (w t) cos (fi t) sin {u t) 



cy{Q) 



< [s^ + cos (11 1) c^] COS (tj t) H sin(fi t) sm{uj t) (? \ z{{S) 



sinwt 



iij 



[l-cos(fit)]scx(0)-sin(f]t)c?/(0) - [s^ + cos (fit) c^] z(0)|, (40) 



where s, c have been defined in ([8]). 

We have used the expressions above to plot the trajectories corresponding to the 
very same initial conditions as in figures of Section |3l The paths are rather different. 

J^..l. Foucault pendulum in 3D 

In Figure Owe represent the exact trajectory corresponding to the Foucault pendulum 
with the same initial conditions as in Figure [31 The left panel presents the xy projection 
and can be compared with the pattern in Figure [3] (left panel) . Notice the difference in 
symmetry. The right panel gives the xz projection. In the two-dimensional case it is 
simply z{t) = 0. 

4-2. Bravais pendulum in 3D 

The Bravais pendulum trajectories are represented in Figure [71 The solid black line 
corresponds to a turn parallel to f2, and the dotted red line is the trajectory with the 
initial tangential velocity reversed. We have used the same condition for the initial 
velocity y(0) as in the planar case. It turns out that after one cycle the oscillation is not 
conical-like any longer. The time delay between both turning modes is also apparent, 
as in the two- dimensions case. 
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Figure 6. Foucault pendulum in three dimensions. Trajectory in the phase plane xy 
(left panel) and in the the phase plane xz (right panel). The solid dot stands for the 
initial position. For the numerical values of the parameters see main text. 
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Figure 7. Bravais pendulum in three dimensions. Trajectories in the phase plane 
xy (left panel) and in the the phase plane xz (right panel). The black solid line 
represents the trajectory when the bob turns in the same direction as 17. The dotted 
line corresponds to a motion started with the reversed initial velocity. The solid dot 
stands for the initial position. For the numerical values of the parameters see main 
text. 



J^.3. General pendulum in 3D 

The more general mode of oscillation of Figure |5] is replicated in Figure |8] for the 
three-dimensional pendulum. The left and right panels give the xy and xz projections 
respectively. As in the two preceding subsections, the differences in the patterns are 
noteworthy. We will discuss about these differences in the final Section. 
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Figure 8. General pendulum in three dimensions. Trajectory in the phase plane xy 
(left panel) and in the the phase plane xz (right panel). The solid dot stands for the 
initial position. For the numerical values of the parameters see main text. 



5. Chevilliet's theorem 

Let us consider the case where we release the pendulum bob at r(0) = (e, 0, A)^, from 
rest: r(0) = 0, with A = / — \/P — e^, and e <^ /. It is then interesting to study the 
trajectory as it is seen from the inertial reference system fixed in space. According to 
( 135|) . the corresponding initial conditions are q{0) = (e, 0, A)^, q{0) = (0, — i7[se+cA], 0)^, 
where s, c have been defined in ([8]). Then (!30|) gives 



q^(t) = e cos{ujt), 

qy(t) = [se + cX] sin(c<;t), 

to 



(41) 



qz(t) = A cos (cut). 

It is straightforward to see that, in three dimensions, the observed trajectory takes place 
on the ellipsoid 



qI + <il 



+ 



? 



^ = 1. 



e2 + A2 (se + cA)2 fi2 
The intersection with the plane g^ = is the ellipse 



Tx 



+ 



Tn 



UJ 



e2 + A2 (se + cA)2 ^2 
whose ratio of axes in the limit e ^ / can be approximated by 



n ( e 

uV^'2l 



(42) 



(43) 



(44) 



since A ^ e^/2/. Thus, the planar projection of the pendulum motion is an ellipse as a 
consequence of the non-vanishing initial (tangential) velocity provided by the rotation 
of the Earth on the bob. Otherwise the trajectory would be on a straight segment. 
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Had we resolved the same initial value problem with the approximated pendulum 
equations in two dimensions of Section |3l we would have found that the trajectory in 
the fixed reference system describes the ellipse 

'£^i.(^Y=l, (45) 



e2 e2 



\sQJ 



whose axes ratio, sQ/u, is determined simply by the ratio of the Earth rotation frequency 
at latitude (3, to that of the natural frequency of the pendulum. This is nowadays 
an almost forgotten result, known in old texts as Chevilliet theorem [HI |12l [T3] . 
Interestingly, we can also recover this theorem by doing a further approximation in 
( H3|) . Indeed, the drastic approximation A = immediately leads to the Chevilliet 
theorem. Thus, the solution in three dimensions (144 p provides a small correction to 
Chevilliet's ellipse f HSj) . We are now ready to accept that the convoluted patterns we 
observe in figures such as Figure |3] is the direct consequence of observing the motion on 
a Chevilliet ellipse from a rotating reference system. 

We can buttress the idea of the geometric origin of the convoluted patterns in the 
phase plane as follows. Consider, in the exact solution in ( 1371) . the subtraction 



S{t) = r{t) - \uj cos ut f (0) + - sin ujt f (0)] . (46) 

We can interpret S as the remaining line of motion that is left once the elliptic piece 
has been removed (i.e, the non-perturbed or inertial part) from the exact trajectory. 
Next, due to the mathematical properties of the cross and scalar products, we get for 
that remainder: Ct ■ S{t) = 0, Vt. Hence, every additional contribution to the trajectory 
originated in Coriolis and centrifugal acceleration terms takes place in a plane orthogonal 
to Cl. There is no correction to the trajectory in the very direction of the Earth angular 
velocity vector Q. 

The flowery nature of the patterns of the pendulum trajectories emerges, therefore, 
from the non-inertial character of the rotating frame. Everything should appear as an 
ellipse in the phase plane if we were able to choose the proper reference system, so to 
speak. 

Finally, neither the solution in two nor the one in three dimensions preserves the 
constancy of the length of the pendulum. This is not a consequence of the non- vanishing 
angular velocity of the reference system but the approximations made to get a linear 
systems of differential equations. It is a price to pay once we have replaced a non-linear 
problem (the physical pendulum) with a linear one (the simple pendulum). Namely, the 
constraint x^ + y'^ + {l — zY = t^ (or, equivalently, qi + <iy + {l — <izY = ^^) which represents 
a sphere of radius /, centered in (0,0,/), is not preserved along a trajectory. Instead 
of a sphere we have the ellipsoid ( H2l) . The same drawback occurs in the approximate 
treatment in two dimensions. There the constraint is x^ + y^ = (or, g^ + g^ = 0), which 
is by no means preserved but replaced with the Chevilliet ellipse constraint (H5|) . Figure 
[9] illustrates this discrepancy in two (left panel) and three (right panel) dimensions as a 
function of time. Both functions are periodic but with different periods. For a pendulum 
length / = 100 m., in the two dimension model the period is 20 sec, and the maximal 
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Figure 9. Foucault pendulum in two (left panel) and three (right panel) dimensions. 
The functions exhibit the non-preservation of the constancy of the pendulum length I 
(preservation corresponds to the value zero in the plots) . 



relative deviation is 1/100. In three dimensions the period is 200 sec, and the the 
maximal relative deviation is 14/100. 



6. Foucault pendulum: Comparison of the 2D and 3D cases 

At this point we have already realized that the trajectories obtained from the 2D- and 
the xy phase plane from the 3D-system of equations do differ. Next we study whether 
these differences are relevant for the analysis of Foucault pendula placed on Earth. 

In his original papers [2] Foucault explains how in Paris the pendulum takes circa 
30 hours to complete one turn, which corresponds to one degree in five minutes. This 
is to say that the museum common visitor observes, at most, that the pendulum bob 
describes a small arc of circumference during her or his visit. In the original setup at the 
Pantheon there was a pile of damp sand disposed in circular shape around the pendulum 
equilibrium point that a stylus attached to the underside of the bob plowed, witnessing 
the accumulative deflection of the trajectory. Modern museums setups rather use falling 
rods or balls instead of sand. 

The solutions ( TT71) and ( IT8l) of the 2D model show that the bob sweeps a circular 
area after a complete turn. This is represented in the left panel of Figure [TOl The black 
dots stand for the location of the bob, sampled every 0.4 sec. (no need to seek any 
especial pattern inside). The pendulum length is 100 m., and Q = 0.000073 rd/sec, is 
the real Earth angular velocity. The bob is released from rest (in the Earth attached 
reference system) at x(0) = 1 m., y{0) = 0. In addition, the right panel in Figure 
[To] represents the corresponding 3D case from (138|) and (!39|) . Here the red solid line 
corresponds to the pattern border in the preceding panel, and it is included for the sake 
of comparison. Definitely, this dual cardioid-like shape is not the observed pattern in 
either museums or reported by Foucault. 
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There is an untold piece in this story that reasonably explains the disagreement. 
Air resistance gradually damps the oscillation of the pendulum. Foucault wrote [2] 
that every five or six hours {i.e., defiection of 60 to 70 degrees) the experiment must be 
stopped in order to restart it from that very precise position with the original oscillation 
amplitude restored. The blue arrow in Figure [10] points out that event and, without 
further information, the receding of the black border could be erroneously ascribed to 
damping. Of course, the air resistance effect is incremental. 

The course of the experiment is then punctuated by restarting and, as a 
consequence, what one sees is the juxtaposition of several copies of the first 60 degrees, 
an (almost) circular sector of the pattern in the right panel of Figure [TOl The overall 
composition mimics the circular pattern. 

Thus damping prevents the experimental observation of the dual cardioid-like 
pattern. The decreasing of the oscillation amplitude is caused by the constructive 
combination of damping and very three dimensional nature of the motion. Indeed, 
Foucault's experiment was not designed to resolve these two contributions but to exhibit 
the defiection of the rotation plane of the pendulum. What we argue is that, in practical 
realizations, both panels in Figure fTOl are not in confiict. 

In 1855 a Foucault pendulum was established as permanent exhibit at Ars et 
Metiers with the occasion of the Universal Exhibition held in Paris. At this time an 
ingenious electromagnetic device was added so that the experiment needed no restarting. 
An electromagnet placed on the ground exerted attraction on the bob only when this 
was approaching the equilibrium point. When it crossed the vertical and the bob was 
receding the action of the electromagnet stopped. A mechanical synchronization device 
restarted the current once the bob was anew in the approaching phase. 

Modern pendula exhibitions incorporate also a driving apparatus. For instance, 
in the Museum of Arts and Sciences in Valencia there is a permanent exhibit. The 
pendulum is driven electromagnetically through a system located in its upper part. The 
device actions are controlled by an electronic automata. The question arises then: to 
what extent is the attractor of the system modified? We will not enter here into this 
subject and will limit ourselves just to point out that the solvable, linear, homogeneous, 
second order system of differential equations with constant coefficients in ([7]) becomes a 
different class of equation: non-constant coefficients, non-homogeneous, or non-linear; 
depending on the type of driving. In general, non-solvable. 

7. Final comments 

The Foucault pendulum as well as its generalizations, are not only beautiful museum 
experiments and excellent tools from a pedagogical point of view (as they contain an 
ideal admixture of mechanics and mathematics), but also they do have far reaching 
applications. The Foucault pendulum-like motion can be viewed as a relative orbital 
motion under a central force, like the relative motion of celestial bodies subject to the 
infiuence of the same gravitational field, or terminal orbital rendezvous and as such the 
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Figure 10. Foucault pendulum. Black dots represent the area swept by the bob 
after one complete turn in the two (left panel) and three (right panel) dimensional 
model. The red solid line corresponds to the border of the pattern in the left panel 
and is given for the sake of comparison. The blue arrow locates 60° deflection of the 
oscillation plane with respect the initial position (large solid black dot). 



subject is still an active topic and has actual relevance. 

The full richness of these systems is sometimes disguised when studied as separate 
entities. The unified treatment we propose here manifestly shows the crucial role 
played by initial conditions in determining the specific features of the system. Here, 
as in spontaneous symmetry breaking, where the Lagrangian possesses a symmetry 
not respected by the vacuum, the equations of motion are identical in every case, the 
selection is done by the initial conditions. 
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